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Abstract 

Recent theoretical advances applied to metamaterials have opened new avenues to design a coat- 
ing that hides objects from electromagnetic radiation and even the sight. Here, we propose a 
new design of cloaking devices that creates perfect invisibility in isotropic media. A combination 
of positive and negative refractive indices, called plus-minus construction, is essential to achieve 
perfect invisibility (i.e., no time delay and total absence of reflection). Contrary to the common 
understanding that between two isotropic materials having different refractive indices the electro- 
magnetic reflection is unavoidable, our method shows that surprisingly the reflection phenomena 
can be completely eliminated. The invented method, different from the classical impedance match- 
ing, may also find electromagnetic applications outside of cloaking devices, wherever distortions 
are present arising from reflections. 
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I. INTRODUCTION 



Rapid progress in electromagnetic and meta-material science research has made it pos- 



sible to envision novel and promising applications in several areas like telecommunications, 
defense, medical imaging and nano-technologies [l-3]. Among them, the design and con- 
struction of cloaking devices that are able to guide electromagnetic fields around any object 
put inside has captured the attention of many researchers. 

In anisotropic media, Pendry et al proposed a coordinate transformation over the elec- 
tromagnetic fields that allows to redirect them at will (4, 5|. The distortion of the fields 
originated by the transformation media is used to generate the new values for e and /i, de- 
fined in anisotropic media. As a result the electromagnetic waves follow the distorted space 
with newly generated e and /i and go around the object and hide it from the sight. 

An alternative approach was introduced by Leonhardt, where a conformal mapping was 
used to design cloaking devices in isotropic media jsl, 7]. The approach describes the prop- 
agation of light based on Hamilton's analogy between the trajectory of rays in media and 
the motion of particles governed by classical mechanics principles. This approach is valid 
when the refractive index profile n{r) does not strongly change over scales comparable with 
the wavelength of light js, 9]. A conformal mapping is then performed between the physical 
space described by a complex field z and an analytical functional w(z) that represents the 



mathematical 
both sheets 



space composed of two Riemann sheets and a branch cut which connected 
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As a result the object remains hidden on the Riemann sheets and 



cannot be detected by an external observer. Although the device does not generate perfect 
invisibility, the distortions can be made relatively small to hide objects that are large in 
comparison with the wavelength. 

Several recent works have extended these ideas to investigate various problems like the 
limitation of the electromagnetic cloak with dispersive material 12J, the extension of the 



bandwidth of electromagnetic cloaks 



space without changing t 



context like acoustics 
carpet cloaking 
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re topology 



13 L the case of transformation media that re-scales 
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15] and have even been applied in a different 



19) . Another interesting application of similar concepts is called 



A carpet mirror can be used to hide an object under it. The carpet 
creates a bump, therefore a distortion can be visible in the reflected image. If a carpet cloak 
is located on top of that bump, it can bend the light so that the distortion is eliminated and 



it gives the impression that the mirror is flat. Researches have extended this concept to the 



3D structure and conducted promising experiments [211 ]. 

In a recent work, we proposed a new design of an invisibility device with isotropic mate- 
rials by combining a conformal mapping and a design composed of alternating positive and 
negative refraction [22|. We call this type of design plus-minus construction. This design 



23 



24j to obtain invisibility effect. However, 



allows us to get around the Nachman theorem 
in order to generate perfect cloaking, we had to fine tune the refractive index profile at the 
branch cut, which represents an inconvenient for a practical realization of the device. In 
contrast, here we show that the plus- minus construction surprisingly cancels the reflection 
phenomena out without any additional modification, even when there is discontinuity of re- 
fractive indices at any boundary, including the branch cut, which usually led to an inevitable 
reflection. 

In this work, we propose a novel method that uses a combination of positive and negative 
refraction indices, called plus-minus construction, that enables us not only to enclose the 
trajectory of light, but also to achieve perfect invisibility in isotropic media, without both 
phase delay and reflection. 

In particular, here we introduce a new construction of cloaking devices by performing 
three successive conformal maps on the space defined by a trivial flat metric. This operation 
results on a physical space R 2 with a non trivial refractive index. A technical advantage of 
using first a flat space metric is that the trajectory of light is trivial in the mathematical space 
and can be easily derived in the physical space through the conformal mappings. Hence, it 
is unnecessary to consider the motion equation for light rays a priori. As we show here, the 
trajectory is readily obtained as the result of operating three successive conformal mappings 
that rotate the space. Interestingly, the proposed dielectric media leads to perfect invisibility 
with absence of both phase delay and reflection. These findings strongly highlight the role of 
the negative refraction in material sciences. Although currently the application of negative 



refraction [25| is confined to the construction of perfect lens [26|, our results show that a 
device consisting of alternating positive and negative refraction (plus-minus construction) 
exhibits extraordinary properties like the perfect cancelation of reflection phenomena. 

Our proposed design has a discontinuity at each boundary between two media with dif- 
ferent refractive indices. It is of common knowledge that discontinuities generate reflections 
leading to distortions. It has been believed that reflections are inevitable in isotropic me- 
dia. Surprisingly, the invented method completely suppresses all reflections that occur at 
the boundary between two isotropic materials having different indices and readily leads to 



perfect invisibility. 

The absolute control of reflection phenomena promises to bring considerable technical 
improvements in diverse areas, encompassing electrical engineering, telecommunications, 
optoelectronics and microelectronics. The practical possibility of guiding electromagnetic 
fields in total absence of reflection phenomena, without using the classical impedance match- 
ing, has relevance to the design and implementation of many optical and electromagnetic 
applications outside of cloaking devices, wherever distortions are present arising from reflec- 
tions. 



II. GEOMETRY AND MEDIA 

Our theory is based on concepts of mapping between non-Euclidean geometry on curved 
space, permittivity e and permeability jx of materials. We briefly review some fundamental 



concepts as follows 



19j: 



In material, the light rays can be bent because of non-trivial permittivity and permeability 
tensors. However, in non-euclidian geometry, light in vacuum can also be bent by non trivial 
metric tensors. We can then design an invisibility device by using a mapping between a 
metric tensor in curved space and the characteristics e and fi defined in physical space given 
a material. The main idea behind this construction is that the trajectory of light in a 
material can be identified by that of non-Euclidean geometry. 

Let e* J and /i 4 - 7 be permittivity and permeability tensors and 7^ be the metric tensor 
in the physical space. And let g^ be the metric tensor in the mathematical space. The 
connection between the metric g^ in the mathematical space and both permittivity e*- 7 and 
permeability ^ tensors in physical space is given by 

a/7 



where g = det(gij) is the determinant of g^-, and 7 = det(jij) is the determinant of 7^. 
Namely, yfg is the volume element of mathematical space, and y/j is the volume element of 
physical space. 



III. ISOTROPIC MEDIA 

In this paper, we focus our attention on isotropic media. 
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Since, in isotropic case, e Jjf and /x Jjf tensors are proportional to the identity matrix, we 
can set 

= r?bij (2) 

where 5ij is Kronecker delta function, and n is an arbitrary scalar function which can be 
interpreted as the refractive index because of the following reason. In isotropic media, the 
metric tensor is given by 

ds 2 = gijdx l dx^ = n 2 (dx 2 + dy 2 + dz 2 ) (3) 
therefore the line element is written as 

ds = n\j (dx 2 + dy 2 + dz 2 ) (4) 
From this, it is clear that the n should be interpreted as a refractive index. 



IV. CONFORMAL MAPPING 

Conformal mappings keep the metric form proportional to the identity as shown in (j2J). 
In this section, we discuss the transformation law under the conformal mappings. Since 
complex analytic function is a conformal mapping, we will consider analytic functions on C 
in our analysis. For z = x + iy, we consider the analytic function as follows: 

z' = f(z) = x'(x,y) +iy'(x,y) (5) 

where i 2 = —1. Then, the Cauchy Riemann equation holds 

dx' dy' 
dx dy 
dx' dy' 
dy dx 

On z' = x' + iy' space, let us set the metric for isotropic media 

ds 2 = n' 2 ((dx') 2 + (dy') 2 ) (8) 
For conformal mapping, the metric is changed as follows. 

ds 2 = n' 2 ((^) 2 + ^) 2 ydx 2 + dy 2 ) (9) 
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(6) 
(7) 



Therefore, the transformation law of refractive index between z space and z' space is given 
by 




We note that this expression can be shown to be the same as the expression in [7| by proper 
variable transformation. 



V. CLOAKING DEVICE CONSTRUCTION AND DESIGN 

We construct our invisibility device by using two ingredients. 

1. First, we operate three conformal mappings which transform the flat space into a 
curved space with non trivial metric. 

2. Next, the negative refraction property of meta- materials which, creates a mirror-like 
image, enables us to enclose the trajectory of light. 



V.l. Conformal map 

We perform three conformal mappings in order to construct our invisibility device. 
Let Z be a complex plane and be identified with physical space. The first conformal 
mapping fx'.Z—> W\ is given by 

wx = h(z) = z + - (11) 
z 

where z e Z, u>i € Wx- Here Wx is a Riemann surface defined by fx- Wx space consists of 
two sheets, the first and second Riemann sheets (Wxf and W\ s ). By the conformal map (TTTj) . 
the outside of the circle (\z\ >= a) is mapped to the first Riemann sheet Wxf, and the inside 
of the circle (\z\ < a) is mapped to the second Riemann sheet Wx s - Hence, Wx = Wxf^W ls - 
The branch cut is given by —2a < Re(wi) < 2a, Im(wx) = 0, which is mapped to \z\ = a in 
Z space. This conformal map is first introduced for invisibility devices in [6]. 

For the second Riemann sheet Wx s , we define the two conformal maps as follows. The 
conformal map f 2 : W ls —> W 2 is 

w 2 = f 2 (wx) = -2a + 4a(^±_^)i (12) 
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where W\ G W\ s and Wi G W 2 . The conformal map fs : W 2 —> W3 is 

^3 = h{w 2 ) = 2a- Aa C W2 A +2<1 f (13) 

4a 

where w 2 G W2 and W3 G W3. The domain and region of map f 2 and fs will be given in the 
next section, in order to make f'2 and fs bijective and restrict both the domain and range 
of f 2 and to the inside the upper part of W\ s , W2 and W3, respectively. Here, we remark 
that (2a, 0) and (—2a, 0) are the fixed points (invariant points) under both mappings f 2 and 
fs- 

On the other hand, we do not perform any transformation on the first Riemann sheet W%f, 
so it is kept invariant. We will discuss later the problem of discontinuity of this mapping on 
the branch cut —2a < Re{w\) < 2a, Im{w\) = between W\f and W± s . 

We summarize the maps fi, f 2 and fs as follows. 

Z 4 Wi = Wi f U W u (14) 

W u ^W 2 ^W 3 (15) 
The inverse function of /1, f 2 and fs are then given by: 

z = f[ 1 (w 1 ) = ± y/wl-4a*) (16) 

U ; 1 = / 2 - 1 K) = -2a + 4a(^±^) 3 (17) 

u* = ^(tos) = 2a - Aa C^ 2 ^ (18) 

We will successively perform three conformal mappings f% , f% , f\ to construct the 
invisibility device. 

We show an illustration of the mapping /3" 1 in Fig. [I] and Fig. |5J The mapping f 2 ~ 1 can 
be visualized in Fig. [3] and Fig. HI 

V.2. Domain and Region of f 2 and f 3 

In this section, we set the domain and region of f 2 and fa, defined in the previous section. 
V.2.1. Domain and region in W 2 

We set the boundary of both the domain of f% and the region of f 2 in W 2 as follows: 

I x = { re ie - 2a|0 < r < 4a, 6 = -tv} (19) 

3 

J 2 = { re ie + 2a|0 < r < 4a, 9 = ^tt} (20) 



Light rays 



(-2,0) 



L 

{2,0) 



FIG. 1. Representation of light trajectory in W3 space, (a = 1) Orange lines are light trajectory 
in W3 space. By the mapping /jj" , the arrow in Fig. [1] rotates 120 degree around the right hand 
side of branch cut (2a, 0), which is shown in Fig. [2] 

The triangle region enclosed by I\, I2 and im(u> 2 ) = is the domain of f$ and the region of 
ji- h and I2 are shown in Fig. [7J In what follows, this region is denoted by /. 

V.2.2. Region in W3 

We set the boundary of the region of f% in W 3 as follows: 



The region enclosed by Ji, J2 and Im(ws) = is the region of fa. J\ and J2 are shown in 
Fig. |6j In what follows, this region is denoted by J. 

V.2.3. Domain in W\ 

We set the boundary of the domain of f 2 in Wi s as follows: 



J2 = h{l2) = {x + iy\2a < x < 6a, y = 0} 



Ji = fs(h) 



(21) 
(22) 



Ki = f 2 (A) = {x + iy\ - 6a < x < -2a, y = 0} 
K 2 = f 2 \h) 



(23) 
(24) 
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FIG. 2. Illustration of light trajectory in W2 space. By the mapping the arrow in Fig. [T] 
rotates 120 degree. Moreover, by the mapping f§ , the light trajectories in Fig. [T] are bent to 
those shown in Fig. [2j 




FIG. 3. Image of light trajectory in W2 space, which is the same as in Fig. [21 but the location 
of the arrow is different, (a = 1) Orange lines are light trajectory in W2 space. By the mapping 
/jf 1 , the arrow in Fig. [3] rotates 120 degree more around the left hand side of branch cut (—2a, 0), 
which is shown in Fig. [H 
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Light rays 




FIG. 4. Representation of light trajectory in W\ space. By the mapping f 2 ~ the arrow in Fig. [3] 
rotates 120 degree more. Moreover, by the mapping /<$" , the light trajectories in Fig. [3] are bent 
more to those shown in Fig. HJ 

The region enclosed by K±, K-i and Im(w\) = is the domain of f2- K\ and K2 are shown 
in Fig. |8j In what follows, this region is denoted by K. 

All orange (continuous), green (dash-dotted) and blue (dotted) lines in Figs. O El HI 
correspond each other by the three conformal maps fx, f2 and f^. 

V.3. Refractive index 

In this section, we will compute the refractive index on physical space Z, by applying the 
transformation law of conformal mappings fflUj) and combining negative refractive index as 
mirror-like image. 

V.3.1. Polar coordinate 

We set the coordinate in each complex space Z, W%, W 2 , W3 space as follows. For each 
space, we set 

z = x + iy (z E Z) (25) 

w x = Xl + i Vl (w 1 eW! = W lf U W u ) (26) 

w 2 = x 2 + iy^ (w-2 E W 2 ) (27) 
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w 3 = x 3 + iy 3 Os e W 3 ) (28) 

We introduce a polar coordinate in W 3 space as follows. 

x 3 = 2a + r 3 cos # 3 (29) 

y 3 = r 3 sin 6 3 (30) 

Here we note that the origin of the polar coordinate is located at the right hand side point 
of the branch cut (2a, 0). 

We introduce a polar coordinate in W 2 space in two ways as follows. 

x 2 = 2a + r 2 cos 9 2 = —2a + r' 2 cos 9' 2 (31) 
2/2 = r 2 sin #2 = r' 2 sin (9 2 (32) 

Here we note that the origin of the polar coordinate (r 2 ,6 2 ) is located at the right hand 
side point of the branch cut (2a, 0). Similarly, the origin of the polar coordinate (r 2 ,9 2 ) is 
located at the left hand side point of the branch cut (—2a, 0). 
Therefore, using ( {TBI , we obtain 

I - <S» 3 < 33 > 

9 3 = 3# 2 (34) 

We introduce a polar coordinate in W\ = Wif U W\ s space as follows. 

xi = —2a + r\ cos 61 (35) 
?/i = ri sin 9± (36) 

Here we note that the origin of the polar coordinate is located at the left hand side point of 
the branch cut (— 2a, 0). 

Therefore, using (jT2l) . we obtain 

1 - (37) 
9 l = 39' 2 (38) 

We introduce a polar coordinate in Z space as follows. 

x = r cos 9 (39) 
y = r sin 9 (40) 

In the next section, we compute the refractive index on the upper area J, I, K in W 3 , 
W 2 and Wis, arid set the refractive index of symmetric area of J, /, K with respect to real 
axis as mirror-like image, respectively, but the sign of refractive index is opposite. 
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V.3.2. Refractive index on W3 



For the region J in the upper area of W3 (the area surrounded by real axis, J\ and J 2 
(See ([2"T]) and fl2"2l ). we set the flat metric as follows 

= (dr 3 ) 2 + (r 3 ) 2 (^ 3 ) 2 (41) 

Considering this, we can set the refractive index in J as follows 

n 3 = 1 (42) 

Therefore in J area (in J C W3), all the light ray trajectories go on the straight line, since 
the metric is flat. (See Fig. E]) 

On the other hand, in the symmetric region of J with respect to real axis in W3, we set 
the refractive index like a mirror like image but the sign of refractive index is opposite 

n 3 = -1 (43) 

Therefore, in the symmetric region of J all the light ray trajectories also go on the straight 
line, since the metric is flat too. However the phase velocity is opposite to n 3 = 1 case (See 
Fig. ED. 



V.3.3. Refractive index on W2 

By applying the transformation law for the conformal map ( fTOl) to / 3 , we can compute 
the refractive index 112 on /, which is the domain of / 3 and the region of fi in W2. In the 
area /, the induced metric is given by 

^ 2 = (n 2 ) 2 {(dr 2 ) 2 + (r 2 ) 2 (^ 2 ) 2 } (44) 

where 

* = (45) 
Therefore in / region, the refractive index is given by 

n 2 = 3(^) 2 (46) 
And in the symmetric region of I with respect to real axis in W2, the refractive index is 

™ 2 = -3(^) 2 (47) 

(See Fig. m>- 
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V.3.4- Refractive index on W\ s 

In the same way, we apply the transformation law for conformal map (fTUj) to f2, we obtain 
the refractive index n\ in K which is the domain of f<± m W± s . In the K, the induced metric 
is given by 

ds 2 = n 2 1 {(dr 1 ) 2 + (r 1 ) 2 (d9 1 ) 2 } (48) 

where 

m = (49, 

Therefore in region, the refractive index is given by 

And in the symmetric region of K with respect to real axis in W\ s , the refractive index is 

- = -<£>"<£>-' ^ 

(See Fig. EJ. 

F.3.5. Refractive index inside the circle \z\ < a in Z 

Finally, we apply a transformation law for conformal map (flOl) to /i. On the region 
R x = f^(K) = {re ie e Z\r < a,-n < 9 < 0} in physical space Z 1 the induced metric is 
given by 

ds 2 = n 2 (dr 2 + r 2 d9 2 ) (52) 

where 

» 2 = (?) 4 (?-)-*(l-^cos20 + 4) (53) 
4a 4a r 4 

Therefore we obtain the refractive index n on the region R\ in physical space Z as follows, 



"=\/(^) 4 (^)-ni-^c O s20+-) (54) 



ro^.ru 4 . 2a 2 „ a 4 , 
^ 4 J. -3 !- cos2# + - 
4a 4a r z 

On the other hand, the refractive index n on the region R 2 = { re ie G Z\r < a, < 9 < n} 
in physical space Z as follows, 



n = _W(!l)4(!l)-!(i_^ C0S 2^+r_) (55) 
" V 4a ; V 4a ; v r 2 r 4; 
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V.3.6. Refractive index outside the circle \z\ > a in Z 



Finally, We set flat metric on the first Riemann sheet of W\ / as follows: 

ds 2 = (dxxf + (dyi) 2 
= (dn) 2 + (ri) 2 (d0O a 



(56) 



Therefore, in all W\f region, the refractive index is given by 



ri\ = 1 



(57) 



Hence, all the light rays go on the straight lines in W\f. 

We apply a transformation law ( ITU]) to f\. Then, for the region R 3 = {re ld G Z\r > a}, 
we obtain 



2a 



a 



n=(l -cos 26 + —) (r>a 



(58) 



Thus, the refractive index n on the region R% in physical space Z is given by 



n 



1 cos 26 + —) (r>a) 



(59) 



V. 3. 7. Summary of Refractive index in Z 



Then, we summarize the refractive index on all area of physical space. 



where 



n = < 



l-^cos2£ + ^) 



(r > a) 



^73(l-^cos2^ + $) (r<a and - vr < < 0) 
1 - ^4 cos26» + 4) (r < a and < < tt) 



(n/4a) 4 / 3 



(60) 



2; = re 

n = \fi{z)+2a\ 

r2 = I/2 o fi(z) - 2a\ 



(61) 
(62) 
(63) 



Therefore, the refractive index in physical space are mixing of positive and negative 
refractive index. We call this type of construction "plus-minus construction" See Fig. [5j 
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FIG. 5. The figure shows the positive and negative regions of refractive index in physical space Z 
(Plus-minus construction). "+" (resp. "— ") means positive (resp. negative) refractive index. 

VI. TRAJECTORY OF LIGHT 

In this section, we compute the trajectory of light. 

The trajectory of light can be obtained by successive application of the three conformal 
maps f$ , from W 3 to Z. In W 3 , all the light rays go along straight lines, since 

the metric of this space is flat. In other words, the refractive index of W3 is 1 or —1 (see 
Fig. |6]). As shown in Fig. [7J in W2, the light rays are bent by the conformal map /3 -1 from 
W3. In wi G Wis, the light rays are bent more by the conformal map f^ 1 from W2 (see 
Fig. [8]). Furthermore, all the light trajectories are enclosed by the negative refraction like a 
mirror-like image. The last trajectory of light is given by applying the conformal map /f 1 
as shown in Fig. [9j 

It is worth noticing that, in this example, we only show the perpendicular incident case 
to branch cut for the matter of convenience. However, we point out that the proposed device 
is invisible for any arbitrary angle of incident light rays. 

On the other hand, the boundary between the device and the invisible region in physical 
space Z can be obtained as follows: 

L=fr\K 2 ) (64) 

Therefore, the area enclosed by L and the symmetric curve of L with respect to the real axis 
is invisible. We can hide any object inside this region as shown in Fig. [9j All the boundary 
represented by red, green, blue lines correspond to each other by the three conformal maps 
U\ f2 l fr 1 in Figs. El El Eland H 
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FIG. 6. Orange (continuous) line represents trajectory of light in W3. Red (dashed) line indicates 
the branch cut. "+" (resp. "— ") means positive (resp. negative) refractive index. All the light go 
along straight lines, since the refractive index on this region is 713 = 1 or n.3 = — 1. Blue (doted) 
line is given by J\ in Eq. ([2T]) , Green (dash-dotted) line is given by J2 in Eq. ([22jh 




FIG. 7. Trajectory of light in W%. Orange (continuous) line is trajectory of light in W2- "+" (resp. 
"— ") means positive (resp. negative) refractive index. Blue (dotted) line is given by 1\ in Eq. f)19|) . 
Green (dash-dotted) line is given by I2 in Eq. (|20p . All the orange, blue and green lines between 
Fig. Hand Fig. [6] correspond to each other by conformal map f^ 1 (see (fToD ) 

VII. PERFECTNESS OF THE INVISIBILITY DEVICES 

The time delay and reflection causes the imperfectness of invisibility devices and many 
other electromagnetic applications. Especially, in our construction, there are several discon- 
tinuous boundaries of refractive indices which may cause reflection (see Eq. ( 160]) ). In spite 
of that, we will show that our device does not generate phase delay and reflection and leads 
to perfect invisibility. 
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FIG. 8. Trajectory of light in W± s . Orange (continuous) line is trajectory of light in W\. "+" 
(resp. "— ") means positive (resp. negative) refractive index. Blue (dotted) line is given by K± in 
Eq. (|19p . Green (dash-dotted) line is given by K2 in Eq. (|20p . All the orange, blue and green lines 
between Fig. [8] and Fig. [7] correspond to each other by conformal map f^ 1 (see (fT7|) ) 




FIG. 9. Trajectory of light in physical space z. Orange (continuous) line is trajectory of light in 
Z. Green (dash-dotted) curve is given by L in Eq. (|64l) . The area enclosed by the green curve 
is invisible. All the orange, green lines between Fig. [9] and Fig. [8] correspond to each other by 
conformal map /f 1 (see PJ) 

VII. 1. No phase delay 

Negative refractive index has opposite sign of phase delay to the positive refractive index. 
The phase delay in the positive refractive index region K in the second Riemann sheet of 
Wis is canceled by the phase shift in the negative refractive index of the mirror like image 
region of K (see Fig. [8]). Therefore, the proposed cloaking device has no time delay because 
of the property of negative refractive index. 
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VII. 2. Boundary reflection 

Discontinuities of refractive index at the boundary between two materials cause the re- 
flection, which may distort the images and precludes the perfect invisibility. In our device, 
discontinuity of refractive index exists in the following two boundaries. 

1. Discontinuity of refractive index at branch cut (see red (dashed) line in Fig. [8]). 

2. Discontinuity of refractive index at the boundary between positive and negative re- 
fractive indices (see blue (dotted) line in Fig. [8]). 

We will show that in both cases, the reflection completely vanishes, using the striking 
features of the plus-minus construction. 



VII. 2.1. Discontinuity of refractive index at branch cut 

Let us consider the boundary between two different material 1 and 2. Let e\ and /ii (resp. 
62 and 112) be permittivity and permeability of material 1 (resp. material 2). Then, the 
refractive index of material 1 (resp. material 2) is given by hi = (resp. h 2 = y/eojl^)- 

By considering the standard electromagnetic theory [23], the reflection and transmission 
ratio are given as follows: 

1. Perpendicular to plane of incidence: 

The transmission ratio from material 1 to material 2 (perpendicular to plane of inci- 
dence) is given by 



Ei f ^cosi + .f^Jl- (?M 2 sin 2 i 



(65) 



where i denotes the incident angle. Similarly, the reflection ratio is given by 

(66) 



P J^-cosi- Ji&Jl- (21)2 sin'' i 
Ei ■ /^cosz+ ,f^Jl- (f 2 ) 2 sin 2 i 



2. Parallel to plane of incidence: 

The transmission ratio from material 1 to material 2 (parallel to plane of incidence) 
is given by 
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Similarly, the reflection ratio is given by 



E r \/ ' i \/ ' 2 V ' "~ ! 



Mi C o S i- JUk Jl- (|l)2 sin 2 i 



Ei / ^cosz + - (?M 2 sin 2 i 



(68) 




If we consider the impedance matching case, 

f-vf (69) 

the previous equations are simplified as follows. For both cases (perpendicular and parallel 
to plane of incidence), we obtain 

E t 2 cos i 



2 cosi 



cosz + Jl - (§^) 2 sin 2 z 



COS % + COS J 

T(i,j) (70) 



where j is outgoing angle of light and 



g r _ cos^-y / l-(fe) 2 sin 2 z 

E * cost + Jl - (|^) 2 sin 2 i 

cos i — cos j 
cos z + cos j 

= R(i,j) (71) 

Let us consider that the light enters from the first Riemann sheet W% / to the second Riemann 
sheet W\s through the branch cut (red (dahed) line in Fig. [H]). In what follows, we call " first 
branch cuf for branch cut from W\f to W\ s and " second branch cuf for branch cut from 
Wis to W\f. At a first look at Eq. (1711) . it seems that there should be some reflections at 
the first branch cut even if we make use of impedance matching. (See red arrows in Fig. 
[TOl ) However, we will show that this reflection is completely canceled out by summing up 
all the reflections that take place at the second branch cut. 
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(2,0) 



FIG. 10. Light rays trajectory in W\ space. Red arrows are reflected light rays at the branch cut. 
Orange and blue arrows show another reflected ray lights. Orange arrows indicate the light rays 
that enter from first Riemann sheet W\ / to the second Riemann sheet W\ s through the first branch 
cut. After that these light rays can also be reflected at the second branch cut. This reflection is 
shown by blue arrows. Blue arrows show that the light rays return to the first Riemann sheet W\t. 
Although figure shows perpendicular incident light rays to the branch cut, the cancelation of the 
reflection is valid for any arbitrary incidence angle. 

Next, let us consider another case that the light enters from Wif to W\ s through the 
first branch cut (orange incident arrow in Fig. [TO]) . This light ray goes around the second 
Riemann sheet (orange arrows in Fig. [TO]) and is reflected at the second branch cut. Then, it 
returns to the first branch cut and goes through the first branch cut again (see blue arrows 
in Fig. [TO]) . In this case, the ratio is given by 



Furthermore, the light rays can be reflected again at the first branch cut and go and back 
two times between the first branch cut and the second branch cut in W\ s , and finally come 
back to Wif through the first branch cut. For the case that the light rays go and back two 
times, the ratio is 



A 1 = T(i,j)R(j,i)T(j,i) 



(72) 



A 2 = T(z,j)i? 3 (j^)m^) 



(73) 



More generally, for going back n times, the ratio is 



Ai = T(i, i)R 2n ~ l (j,i)T(j,i) 



(74) 
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The total sum of all contribution is 

Atat = T(i,j){R(j,i) + R 3 (j,z) + R 5 (j,i) + . . .}T(j,i) 

oo 

= T^j){J2R 2n -\j^)}T(j^) (75) 
i=i 

We insert Eq. f JTOj) and Eq. ( ITT]) into the previous equation, and the total sum of the 
contributions can be computed as follows. 

cos j — cos i . . 

= - 76 

COS I + COS J 

Unexpectedly, the sum of all contributions cancel the original reflection R(i,j) of Eq. 
( I7T1) . Therefore, there is no reflection at the branch cut at all. The perfect cancelation of 
reflection phenomena is surprising itself and can be seen as a consequence of the proposed 
plus-minus construction. 



VII. 2. 2. Discontinuity of refractive index at the boundary between positive and negative refrac- 
tion indices 

Next, we discuss about the discontinuity of the boundary between positive and negative 
refraction indices. (See blue (dotted) line in Fig. [8]) 

If we set ei = — e 2 and Hi = — /i2 at the boundary, we obtain 

hi = -n 2 (77) 

In this case, by considering Eq. fl66|) and Eq. fl68|) (or Eq. flTTj) ). we see that there is no 
reflection at all for all incident angle and polarization at the boundary between positive and 
negative refractive materials. 

f = ° (78) 
As a conclusion, there is no reflection for both the boundary at branch cut (red (dashed) 
line in Fig. [8]) and the boundaries between positive and negative refractive materials (blue 
(dotted) line in Fig. [8]). Therefore, this construction is perfect in theory, with no time delay 
and complete absence of reflection. 

It is worth noticing that the perfect invisibility effect (i.e., no time delay and no re- 
flection) is independent of the precise refractive index profile. In order to create a perfect 
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cloaking, it is enough to enclose the light trajectory in W\ s by a proper symmetric positive 
and negative refractive index (plus- minus construction). Therefore, the current proposed 
construction is just one example because many other plus- minus constructions that lead to 
perfect invisibility are possible using the same principles. 



VIII. SUMMARY 

We have presented a new design of perfect cloaking devices. The design involves a new 
method as well as a construction based on negative refraction material. The novel method- 
ology shows that we do not need to use the equation of motion for light rays. The trajectory 
of light is obtained easily by conformal mapping rotation, which simplifies the theoretical 
construction of the model. Second, the negative refraction material is crucial to achieve per- 
fect invisibility. The results not only indicate that there is no time delay but also show that 
the reflection phenomena is completely canceled out thanks to the plus-minus construction. 

The finding that reflection is completely suppressed suggests that the plus-minus con- 
struction may be potentially applied to various areas and technologies outside of purely 
invisibility oriented devices. In our view the fundamental roots of the complete cancelation 
of light reflection are worthy of further exploration. The absolute control of reflection phe- 
nomena promises to bring significant technical improvements in diverse areas, encompassing 
electrical engineering, telecommunications, optoelectronics and microelectronics. The prac- 
tical possibility of guiding electromagnetic fields in absolute absence of reflection phenom- 
ena has relevance to the design and implementation of many emerging applications such 
as lenses, medical imaging, microwave passive devices, radar/defense systems and wireless 
communications, from cell phones to Global Positioning Systems. 

Impedance matching is the traditional method to reduce the reflection for various de- 
vices. Here, we have demonstrated that a design based on the plus-minus construction can 
completely eliminate the reflection that occurs at boundaries between two materials hav- 
ing different indices of refraction, even though both indices are positive. Taken together, 
the invented method suggests a new way of reducing the reflection, beyond the classical 
impedance matching, that may lead to improve the efficiency of a large variety of optical 
and electromagnetic devices. 
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